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Abstract — In this paper, we investigate joint optimal relay 
selection and resource allocation under bandwidth exchange (BE) 
enabled incentivized cooperative forwarding in wireless networks. 
We consider an autonomous network where TV nodes transmit 
data in the uplink to an access point (AP) / base station (BS). 
We consider the scenario where each node gets an initial amount 
(equal, optimal based on direct path or arbitrary) of bandwidth, 
and uses this bandwidth as a flexible incentive for two hop 
relaying. We focus on a -fair network utility maximization (NUM) 
and outage reduction in this environment. Our contribution is 
two-fold. First, we propose an incentivized forwarding based 
resource allocation algorithm which maximizes the global utility 
while preserving the initial utility of each cooperative node. 
Second, defining the link weight of each relay pair as the utility 
gain due to cooperation (over noncooperation), we show that the 
optimal relay selection in a-fair NUM reduces to the maximum 
weighted matching (MWM) problem in a non-bipartite graph. 
Numerical results show that the proposed algorithms provide 20- 
25% gain in spectral efficiency and 90-98% reduction in outage 
probability. 

I. Introduction 

The benefits of cooperative communications [ 1 ] have led to 
significant research in relaying and forwarding [2 |. However, 
forwarding always incurs some costs, e.g., power and/or delay. 
There have been few works that focus on the explicit cost 
of forwarding. Existing cooperative communications literature 
include several incentive based mechanisms to encourage the 
forwarder nodes for cooperation. These techniques include 
pricing (3), reputation |4| and credit |5| based cooperative 
forwarding. However, these mechanisms require a stable econ- 
omy or a shared understanding of what things are worth and 
become unrealizable in dynamic wireless networks. 

In light of this, the authors of (6] developed a bandwidth 
exchange (BE) enabled incentive mechanism where nodes 
offer a portion of their allocated bandwidths to other nodes 
as immediate incentives for relaying. They used a Nash 
bargaining solution based resource allocation and a heuristic 
relay selection policy in their work. In this work, we focus 
on the distributed joint optimal relay selection and resource 
allocation in the a-fair NUM and outage probability reduction 
of a BE enabled network. 

We consider an N node autonomous network where each 
node receives an initial amount (equal, optimal based on direct 
path transmission or arbitrary) of bandwidth and connects 
directly to the access point (AP) / base station (BS). We 
consider a frequency division multiple access system where 
all nodes transmit at the same time with different bandwidth 



slots. In this context, we focus on a two-hop incentivized 
cooperative forwarding scheme where a sender node provides 
bandwidth as an incentive to a forwarder node for relaying its 
data to the AP/BS. We first prove the concavity of the resource 
allocation problem and then show that the optimal relay 
selection problem in a-fair NUM reduces to the classical non- 
bipartite maximum weighted matching (MWM) algorithm |7J . 
Using the distributed local greedy MWM [8|, we propose a 
simple distributed BE enabled incentivized forwarding pro- 
tocol. We also show that the outage probability reduction 
problem reduces to the bipartite maximum matching algorithm 
in this context. Numerical simulations show that the proposed 
algorithm provides 20-25% spectrum efficiency gain and 90- 
98% outage probability reduction in a 20 node network. 

A. Related Work & Our Contributions 

Our contributions in this paper can be summarized as 
follows. First, we consider incentivized relaying in a network 
where each node has been allocated an initial amount of 
bandwidth. Previously, the authors of (9) and flO) considered 
incentivized forwarding in a cognitive radio network where 
only the primary users initially receive resources and later 
transfer some of their resources to the secondary users as 
incentives for relaying. In contrast, our work focuses on 
distributed incentivized two-hop relaying in an autonomous 
network where a centralized algorithm might be infeasible due 
to the associated long estimation delay and high complexity. 

Second, our proposed decode & forward (DF) BE enabled 
resource allocation maximizes the summation of the utilities 
while preserving the initial utilities of the individual nodes. 



Previously, the authors of 1 1 1 1 considered BE from a simpler 
two hop relaying perspective. The authors of [12] proposed 
a similar half duplex DF relaying approach. However, they 
considered a commercial relay network where the relay did 
not have its own data fT2) . To the best of our knowledge, our 
proposed BE based resource allocation algorithm has not been 
investigated before. 

Third, our definition of link weight in the use of MWM is 
different from that of existing literature. Inspired by the semi- 
nal work on maximum weighted scheduling [13], most of the 
works on MWM based scheduling have defined link weights 
as the differential backlog size of that particular link 1 14], 1 15 1. 
Based on this definition and using the network layer capacity 
perspective, the MWM algorithm of these works finds the set 
of links that will be activated at each slot fl6) , (17). However, 




we adopt an information theoretic capacity perspective in this 
work. We define the link weights of the MWM graph as the 
utility gain that a DF relaying enabled cooperative pair offers 
to the system. Therefore, the 'matched' nodes of the MWM 
algorithm communicate to the AP using a DF cooperation 
strategy whereas, the 'unmatched' nodes transmit to the AP 
without cooperating with any other node. In this regards, 



our relay selection approach is closer to the work of [18 1 
where the authors defined link weights of each cooperative 
pair as the energy savings of cooperation over noncooperation. 
However, [ 18 1 considered energy minimization from a bit error 
rate (BER) perspective whereas, we focus on a-fair NUM from 
capacity perspective. 

This paper is organized as follows. Section [II] and III 
show the proposed incentivized system model and the de- 



sign objective respectively. Section IV describes the problem 
formulations and solutions of incentivized forwarding based 
resource allocation and relay selection. After providing the 
simulations results in Section [V] we conclude the work in 
Section El 

II. System Model 

We consider the uplink of an A node single cell FDMA 
network. Let V = {1, 2, • • • , A} denote the set of A nodes 
that transmit data to the BS (node 0). Each node uses the total 
time slot. Node i £ Vis initially allotted a bandwidth of W\ n . 
Let pij and Rij denote the link gain and achievable rate of 
the ij path respectively. For a given node i G [1, A] with W\ n 
bandwidth and direct link gain pio, the achievable throughput 
is: 



R<r=Ri = wriog 2 1 



wr 



(i) 



Here, Pi, ma x denotes the maximum transmission power of 
node i and R\ n is expressed in bit per second (bps). 

In BE, nodes perform two hop half duplex DF cooperative 
relaying. The forwarder node hears the sender nodes' data and 
relays some of that data along with transmitting its own data 
to the AP. Since each node is a power constrained device, the 
forwarder nodes' data rate may drop if it continues to transmit 
with the same bandwidth. Therefore, the sender node delegates 
some of its bandwidth to the forwarder node as incentives for 
relaying. We consider one forwarder for one sender and vice 
versa to reduce the relay searching complexity. Let SJ- = 



TABLE I 
Summary of used notations 



Notation 


Meaning 


N 


Total number of users 


P 


Maximum transmission power 


pij 


Gain of the ij link 


w m 


Initial bandwidth of node i 


w be 


Node i's bandwidth in BE 


R T 


Initial rate of node i 


R? 


Node j's rate in BE 


Rij 


Achievable rate in the ij link 


V 


Set of all nodes 


V 


Set of nodes that transmit without cooperation 


ST 


Set of sender-forwarder pairs 


SFi 


i th sender-forwarder pair (s;,/i) 



{STi,--- ,ST K } = {(si,/i),(s2,/ 2 ),--- ,( s kJk)} de- 
note the sender-forwarder pair set, i.e., fi relays s/s data 
along with transmitting /j's own data. Let T> = di, (fe, • • ■ , di, 
denote the direct set, i.e., the set of remaining nodes that 
transmit data without cooperation. Note that, K and L are 
variables and further, 2 * K + L = N. We assume pairwise 
bandwidth constraint in this work. 

The left and right sides of Fig. [T] show the considered 
direct transmission model and the proposed BE model. In BE, 
node 2 relays data for node 3. Node 3 delegates a amount 
of bandwidth to node 2 as incentive for relaying. Let W\ e 
represent the bandwidth of node i in the BE scenario. Now, 
Wt = Wi n + a and I^| e = W\ n - a. Node 4 and 5 operate 
in the same manner. 



Nodes don't do power allocation among different streams 
in our framework. Since capacity is a non-decreasing function 
of transmission power, each node utilizes the maximum trans- 
mission power in its allotted bandwidth slot. Without loss of 
generality, we assume Pi_ max = PVi € V in the subsequent 
analysis. 

1) Rate analysis in the BE scenario: : Let R\ e represent 
the achievable rate of node i in the BE scenario. The right side 
of Fig. [2] shows the interaction between a sender node s and 
a forwarder node /. If node s transmits to node / and the BS 
separately using W^ e bandwidth, the achievable throughput in 



the respective paths are: 

R s f = 
R s o — 



W be loKofl 



If / transmits to the BS using Wi e bandwidth, 



R f0 = W be \og 2 (l 



(2) 
(3) 

(4) 



Assuming p s t > p s0 , it is easily seen that R s f > R s0 . 
Due to the nature of the wireless environment, when sender 
s transmits, both / and BS hear it. If node s transmits 
at rate R s f, node / can decode it properly. The BS also 
receives the same signal but can't decode it properly since 
Rsf > Rso- However, node / can forward R c bits to the BS 
to resolve the BS's uncertainty about node s's data. Node / 
also transmits its own data, R be , to the BS. The information 
theoretic generalization of the maximum-flow-minimum-cut- 
theorem 1 19 1 provides the following relationship between these 
achievable rates, 



R, 



< 
< 



mm(i? s /, R s0 + R c ) 



R 



fo 



(5) 
(6) 



The codebook design procedure to achieve these rates is 
summarized in Appendix [A] A detailed description can be 
found in (19), (20). 

Using this rate analysis, we focus on the distributed joint 
optimal bandwidth allocation and relay selection in maximiz- 
ing the summation of the a-fair utilities of the network. 

a-fair utility: 

The a-fair utility is defined for any a G [0, oo), as (21), 



(R^L if 
U a {R) = I l ~J \ 

\log(R), if 



a ^ 1 
a = 1 



(7) 



where R represents the rate of the user. Summation of a-fair 
utility functions takes the form of different well known utility 
functions, e.g., sum rate maximization (a = 0), proportional 
fairness (a — 1) and minimum rate maximization (a = oo). 

III. Design Objective 

We begin with a focus on a-fair NUM of the overall network 
through optimal bandwidth and rate allocation for all possible 
sender-forwarder pairing sets. 



Problem I 



max . 



E 



U a (R b d ' 



be \ 



E ( 



U a (R b f) + U a {R be )) (8a) 



Wf + W be < W} n + W™ V (s, /) G SF (8d) 



Wf, W° e >0V(s,/) G ST 



(8e) 



V C V , SF G V x V , SFi n ST 3 =$Vi^j (8f) 



STinV = 0Vi G \1,K] 



SJ 7 ! U SF 2 • • • U SF K U V = V 



(8g) 



(8h) 



Variables V, ST, R b s e , Rf,W be , W 



■ f 



Here, R b d e = R d n V G D. Therefore, the rates of the direct 



node set are not optimization variables. Equation ( [85) denotes 
that the rate of the sender and forwarder lie in the convex hull 
of the allocated bandwidth. The details of this convex hull 
has already been explained in the system model. It will also 
be mentioned in the next section. Eq. (8c| > represents that the 
sender and the forwarders' rate cannot drop below their initial 
rates. Equation (Sd} shows that the total bandwidth used by the 
cooperative pair is constrained by the summation of the initial 
bandwidths allocated to the individual nodes. Equation (|8f|i- 
( |8"h"l ) represent the relay selection constraints. Equation (8Tj ) 
shows that the direct and sender-forwarder sets are all subsets 



of the overall set. Eq. ( 8g i denotes that the sender-forwarder 
pairs and direct set cannot have any common nodes. Eq. ([8KJ> 
represents that the union of the pairs and the direct set form 
the overall set V. 

The solution of the above optimization problem depends 
on the selected sender-forwarder and direct node set and 
the corresponding bandwidth and rate allocations. Hence, it 
involves an exponential number of variables and constraints. 
In the rest of the paper, we focus on solving this problem. 

IV. Optimization Problem Solution 

A. Modified Optimization Problem 

Let U a {R to t) = J2 ie v U a (Rf l ) denote the summation of 
the initial utilities of the nodes. For a fixed SJ-, U a (Rtot) can 
be expressed in the following form: 

U a (R to t) 

Yu a {RT) 



E ua ^) 



(U a (Rf) + U a (Rl n ))(9) 
(U a (R l J l ) + U a (Rl n ))lO) 



s.t. (R b t e ,R\ e ) G conv(W b f e ,W be ) V(a, /) G SF (8b) Equation m follows from eq. jHEJ. Equation {TO} uses the 

1 1 foot that j?oe _ r/invjj ^ n c„l, tr ,„ t ;„„ Tja(U \ frr.™ 



Rf > Rf 



RT >i?rV(s,/) G ST 



(8c) 



fact that R b d e = R d n \J d G D. Subtracting U a {R tot ) from 
the objective function of I, we find the following optimization 
problem: 



Problem II 

(U a {Rf) + U a (R b s e )-U a (K l f n ) + U a {RT)) (11a) 

s.t. (R b f e ,Rl e ) G conv(W b f e ,W be ) V(s,/) G 5J 7 (lib) 

i?^ e > Rf , R b s e > R™ V (s, /) G ST (1 lc) 

W) e + W s be < W} n + WT V (s, /) G 5J- (1 Id) 

W be , W be > V (s, /) G ST (lie) 

X> c V , ST G V x V , n = Vi^ j (llf) 

5J 7 < n2? = 0V< G (llg) 



J2$? > Rf , R be > Rf V(s,/) e 5J" (12c) 

W) e + W s 6e < W} n + W l s n , W b f e , W be > V(s,/) G 5J" 

(12d) 

Variables R b s e , R b f e , W be , Wf (12e) 

Now, due to the pairwise bandwidth constraint of ( |12d| i, the 
bandwidth allocation in one cooperative pair does not affect 
other nodes. Therefore, problem III is just the summation 
of K independent three node (sender, forwarder and BS) 
resource allocation problems and can be decomposed into the 
subproblems. Hence, we now focus on an arbitrary sender- 
forwarder pair (s, /) and describe the resource allocation 
problem formulation in this pair. 

Problem IV 

max. U a {Rf ) + U a (R be ) - U a {Rf) - U a {R\ n ) (13a) 



ST!UST 2 -- - USTk^V = V (llh) 

Variables V, ST, R be , R b f, W be , Wf 

The inclusion of constant terms in the objective function 
does not change the optimal variables of an optimization 
problem [22 1. As a result, the same set of sender-forwarder 
pairs maximize both problem I and II. We will focus on solving 
problem II in the subsequent analysis. The optimal variables of 
problem II will directly lead to the optimal solution of problem 
I. 

Problem II depends on both relay selection and resource 
allocation. The very nature of the design objective allows us 
to split the optimization formulation into the following two 
parts: 

• For any fixed set of sender-forwarder pairs, perform DF 
based optimal rate and bandwidth allocation. 

• Choose the relay set that maximizes the summation of 
the a-fair utility of the nodes. 

We, at first, focus on resource allocation in a fixed sender- 
forwarder pair set ST and direct node set T>. Later, we will 
show the optimal sender-forwarder pair selection policy. 

B. Optimal bandwidth and rate allocation for a fixed sender- 
forwarder set 

The optimal resource allocation problem for fixed ST and 
T> takes the following form: 

Problem III 

J2 (U a {R b f) + U a (R be )-U a (Ry)-U a (Rl n )) (12a) 
s.t.(R b f e ,R b s e ) € conv(W b f e ,W be ) V (s, /) G ST (12b) 



s.t. i? s/ <W^log 2 (l + ^) 
^o<^ e log 2 (l + ^°) 
R c + Rf<W b f e\o g2 (l + ^^) (13b) 

R be < mm(R sf ,R s0 + R c ) , R b f > Rf , R be > R™ (13c) 

W b / + W be < W} n + Wl n , Wf, W be > V(s,/) G ST 

(13d) 

Variables R be , R be , W be , Wf , R c 

Equation ( |13b| i and ( | 1 3c[ > show the convex hull of the allotted 
bandwidths and the achievable rates. They have also already 
been described in the system model. 

Lemmma 1: Problem IV is a concave maximization prob- 
lem. 

Proof: U a {Rf) and U a {R™) are the utilites of the initial 
data rates and constants, in terms of the optimization variables. 
The concavity of a-fair utility functions and the capacity 
expressions can be easily shown [22]. The minimum of linear 
(concave) functions is concave. Thus, the objective function is 
concave and the constraints are convex or linear in terms of the 
optimization variables. This proves the concavity of Problem 
IV. ■ 

Problem IV can be solved using standard convex optimiza- 
tion algorithms, e.g., interior point methods p2) . 

Lemma 2: {R be , R b f , W be ,Wf,R c ) 

(Rl n ,Rf,W* n ,W} n ,0) is a feasible set of variables of 
Problem IV. 



Proof: Let Wf = Wf 1 , W b s e = W™, R c = 0. Then, 
R sf = Wriog 2 (l + ^l) 
R sO = W ™log 2 {l + ?^0-)=R? (14) 

VV S 

Rf = mm(R sf , R s0 + R c ) = min(i? s/ , R s0 ) = R™ 

This is a feasible solution of problem IV. ■ 

Thus, if Node s and / use their initial bandwidths and if 
the forwarder / does not relay any data of the sender s, both 
sender and forwarder will continue to transmit at their initial 
rates. The corresponding feasible solution for these variables 
can be found as follows, 

U a (Rf) + U a (R b s e ) - U a {Rf) - U a (R l s n ) 

= U a {Rf) + U a {RT) - U a (Rf) ~ U a (R l s n ) = 

Since, problem IV is a concave maximization probelm, the 
optimal solution will not drop below J22) . Hence, the 
proposed BE enabled relaying scheme will perform at least as 
good as the initial allocation and seek to maximize the global 
utility while preserving the initial rates of the individual nodes. 

Lemma 3: Problem III is a concave maximization problem. 

Proof: Problem III is the combination of K disjoint concave 
maximization problems and hence it's concave. ■ 

C. Optimal sender-forwarder set selection 

Let g*j denote the solution of problem IV, i.e., the optimal 
gain obtained through cooperation of (s, /). Let, sjr 
represent the optimal solution of problem EH, i.e., the optimal 
cooperation gain for the selected sender-forwarder pair set 
SJ-. Then, C/p SJ r = Y^( s ^esrdsf- Therefore, the sender- 
forwarder pairing set selection part of problem II is equivalent 
to finding the set of pairs that maximize the gain in utility 
through cooperation, over noncooperation. It can be written in 
the following form: 

Problem V 

argmax ^ g* f (15a) 

5JeVxV, SFi n ST 3 = ®Vi^j, (15b) 

Now, consider a graph Q — (V,£) where the vertices V 
represent the set of N nodes under consideration and £ denote 
the edges between these nodes. Define the edge weight of any 
pah- by U a {R he ) + U a (Rf)-U a {R l j n )-U a {R™), i.e., 
the difference, in terms of utility, between the cooperation and 
non-cooperation scenario. The optimal sender-forwarder pair- 
ing set selection of problem V is equivalent to finding the set 
of pairs that maximize the difference between cooperation and 
noncooperation utility. Hence, the optimal sender-forwarder 
pair selection problem can be reduced to the problem of 



finding the set of pairs that maximizes the link weights 
mentioned above. Thus, the optimal relay selection converges 
to the classical nonbipartite MWM problem. The nonbipartite 
MWM algorithm has been summarized in Appendix [B] A 
detailed description can be found in (7). 

Fig. [3] illustrates the application of MWM in the sender- 
forwarder pairing selection. The left figure of Fig. [3] denote 
the initial scenario where node a, b and c transmit through the 
direct path and transmit 6, 3 and 7 bits respectively. The three 
figures in the middle show the obtained rates for different 
sender-forwarder pair selection. The second figure from the 
left shows that b and c transmit 4 and 9 bits respectively 
through BE enabled DF cooperation. Thus, the utility gain of 
cooperation, over noncooperation, is 3 bits. The middle figure 
and the 2nd figure from the right represent the cooperation 
scenarios of (a, b) and (a, c) respectively. The rightmost figure 
represents the edge weights of each cooperative pair in terms 
of the utility gain of cooperation, over noncooperation. The 
MWM algorithm will select node (b, c) as the cooperative pair 
and node a will transmit without cooperation. 

Centralized nonbipartite MWM can be solved optimally 
in 0(N 3 ) time |7J. Our proposed distributed incentivized 
forwarding is based on the distributed local greedy MWM 1 8 ] 
(0(N 2 )) and is described below. The distributed MWM 
algorithm [ 8 1 is summarized in Appendix |C] 

D. Distributed BE incentivized forwarding protocol 

• Focus on an arbitrary node, node v. Node v sends training 
symbols to the AP and obtains its own direct channel, 
p v o, through feedback. Node v is initially allocated W™ 
bandwidth and transmits at R™ rate. 

• Let node a be a neighbouring node of v. Due to the 
nature of wireless channels, node u receives node v's 
channel estimation training symbols and finds the inter- 
node channel gain, p uv . 

• v sends an omnidirectional signal containing p v0 , W™ 
and R™ to the neighbouring nodes. 

• Node u may relay node v's data if min(p uv , p u o) > p v o. 
Thus, v knows its potential forwarders or senders. 

• v solves problem IV for the suitable neighbours. Thus, 
each node knows its adjacent link weights. 

• v solves the distributed local greedy MWM algorithm 
of (8). 

• The 'matched pairs' allocate resources among them- 
selves. The 'unmatched' nodes transmit without coop- 
eration. 

E. Outage probability reduction in BE 

We define outage probability as the ratio of the number of 
nodes that do not get minimum data rate to the total number 
of nodes. We assume that each node in the network starts with 
an initial amount of resource. Depending on the resource and 
link gains, nodes fall in the following two groups: 

• Outage group: Node that cannot meet the minimum 
required rate with initially allocated resources. 






b b b 

Fig. 3. MWM in BE Enabled Relay Network 



• Non-Outage group: Node that can meet the minimum rate 
with initially allocated resources. 

The outage probability reduction problem can be defined as 
providing minimum data rate to the most number of users in 
the outage group, while maintaining the minimum data rate 
of the nodes, in the non-outage group. MWM based matching 
and pairwise resource allocation based incentivized two-hop 
forwarding can help in this case. We propose the following 
scheme in this regard: 

• Each node in the outage group solves the pairwise sum- 
rate maximization, with minimum rate constraints, for 
each of its neighbouring node. Nodes can solve sumrate 
maximization by plugging a = in the a-fair utility 
function. 

• If the node in outage can maintain minimum data rate by 
pairing with the forwarder, i.e., the non-outage node, we 
assume that an edge exists between these nodes. 

• The relay selection problem in outage probability reduc- 
tion becomes maximizing the number of edges in the 
network. This reduces to the maximum matching (MM) 
algorithm in a bipartite graph J23) . 

Our focus here is to maximize the number of users that 
receive minimum data rate, not to maximize any utility func- 
tion. That's why, we use MM, instead of MWM, in this part 
of the work. Besides, cooperation between two nodes in the 
non-outage group do not change the outage probability of the 
network. Therefore, we consider a bipartite graph by dividing 
the graph into outage and non-outage group. 

V. Numerical Simulations 

We assume equal initial bandwidth allocation in all of our 
simulations, i.e., nodes start with equal bandwidth. However, 
our work is readily applicable to the scenario where nodes 
start with optimal bandwidth allocation (based on direct path 
transmission) and then use bandwidth as incentives for two 
hop relaying. 

Fig. [4] and Fig. [5] compare the performance of BE re- 
laying with that of direct path transmission in the sumrate 
maximization and minimum rate maximization of a 3 node 
network (sender, forwarder and BS). Both sender and for- 
warder initially receive 10 MHz bandwidth and transmit uplink 
data to the BS. The forwarder node is placed in the straight 
line that connects the BS and the sender node. The distance 
between the forwarder node and the BS is kept fixed at 150m, 
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Fig. 4. Sumrate Maximization in a 3 Node Network, 10 MHz per Node, P 
= 100 mW, Near node-AP distance = 150m 



Forwarder rate (direct) 
Sender rate (direct) 
Forwarder rate (BE) 
Sender rate (BE) 
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Fig. 5. Minimum Rate Maximization in a 3 Node Network, 10 MHz per 
Node, P = 100 mW, Near node-AP distance = 150m 



whereas, the distance between the BS & the sender node is 
varied. In these two simulations, we assumed the link gains 
to take the form, = kd~ 3 where is the distance 
between the i th and j th node, k is the proportionality constant 
that also captures the noise spectral density and is set to 
k = 6 x 10 6 MHz * m 3 /mW @. 

The sumrate maximization objective based plot of Fig. [4] 
shows that BE relaying improves the rate of the forwarder 
(near user) while ensuring that the sender's (far user) rate 
does not drop below its initial value. This increase in the 
forwarders' rate gets reflected in the sumrate gain of the 
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- BE (Distributed) 
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tributed algorithms in Fig. [6] and Fig. [7] In the simulation of the 
distributed algorithm, we assumed that each node can only talk 
to its neighbours that are located within 500m from the sender 
node. Fig. [6] shows that the centralized and decentralized 
algorithm improves the spectral efficiency by 25% and 20% 
respectively. The performance of the distributed algorithm will 
improve if we allow each node to talk to neighbours with 
greater distances. 

Fig. [7J shows that BE enabled relaying provides cooperative 
diversity and significantly reduces the outage probability (90- 
98%). Thus, BE can be used to extend the coverage in an 
autonomous network. 



Fig. 6. Spectrum Efficiency in an N Node Network, 1 MHz per Node, P 
= 20 dBm 
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Fig. 7. Outage Probability in an N Node Network, 1 MHz per Node, P = 
20 dBm 

network. On the other hand, the minimum rate maximization 
objective based plot of Fig. [5]shows that BE relaying improves 
the senders' (far user) rate while ensuring that the forwarders' 
(near user) rate does not drop. This diverse contribution 
of BE relaying comes from the problem objectives of the 
respective simulations; sumrate maximization (a = 0) is the 
most efficient allocation whereas minimum rate maximization 
(a — oo) is the most fair one. Therefore, the use of any 
a 6 (0, oo) would have increased both users' rates in this 
simulation scenario. 

Fig. [6] and Fig. [7J show the performance of BE relaying in 
an TV node network. In these simulations, we assumed that 
links are under independent Rayleigh fading and the link gain 
in each slot is an independent realization of Rayleigh random 
variable. This means that the link gain, p^, is exponentially 
distributed, p(pij) = -^exp^—j^) where pij — fed" 3 and 
k = 6 x l0 6 MHz * m 3 /mW. We consider a circular cell 
of 800m radius. The AP is located at the center, whereas, 
the nodes are placed randomly in the cell. We considered 
transmission scheme much like the one used in mobile Wimax. 
Each node is preassigned 20 dBm transmit power and 1 MHz 
bandwidth. We used the matching code of |24| to implement 
the MWM algorithm in Matlab. 

We showed the performance of both centralized and dis- 



VI. Discussion 

In this paper, we considered joint optimal relay selection and 
resource allocation in the a-fair NUM and outage probability 
reduction of a BE network. Our proposed resource allocation 
formulation maximizes the global utility of the cooperative 
pair while preserving the initial utilities of each individual 
node. We showed that the relay selection part of the a-fair 
NUM problem reduces to the nonbipartite matching algorithm. 
Numerical simulations suggest that the proposed BE enabled 
relaying provides 20-25% spectrum efficiency gain and 90- 
98% outage probability reduction in a 20 node network. 

Our work can be really advantageous in the following 
scenario: the nodes start with a direct path based centrally 
allocated optimal resources. Thereafter, nodes employ the pro- 
posed algorithm to improve the system performance through 
distributed relay selection and pairwise resource allocation. 
This scheme saves a huge amount of signalling overhead, an 
inherent drawback of centralized optimal two hop forwarding, 
at the cost of lower performance. 

The proposed algorithm considers one forwarder for one 
sender and vice versa. The generalization of this algorithm 
to the multiple sender-forwarder scenario is an area of future 
research. We also plan to focus on discrete subcarrier exchange 
in our future works. 
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Appendix A 
Codebook Design in Proposed DF Relaying 

Let's consider two codebooks VV and C that consist 
of 2 Rs f and 2 Rc codewords respectively. Assume, VV = 
{wx,W2, ■ ■ ■ ,w 2 R sf } and C = {c\, C2, • • • , c 2 k c }. Here 
Rsf > Rc- Consider a partition S = {81,82, ■•• ,S 2 Ro} 
of VV, i.e., VV has been partitioned into 2 Rc cells. Each cell 
Si contains 2 R "f~ Rc codewords of VV. Assume a one-to-one 
correspondence between C and S, i.e., each codeword of C 
represents one particular cell of S. 



The BS (node 0), sender s and forwarder / get the code- 
books off-line. At the beginning of transmission, sender s 
sends a codeword Wi from W using R s f bits. The forwarder 
node decodes the codeword correctly. However, since R s f > 
R s q, the BS cannot decode it correctly. The BS has a list 
of possible codewords of size 2 Rs f~ Rs0 . Now, the forwarder 
/ finds the cell Si where w, lies and sends using R c 
bits. The BS receives Cj and intersects Si with the list of 
possible codewords. If R c > R s f — R s q and R c < Rfo, 
this half duplex DF cooperation completely removes the BS's 
uncertainty about Wi fl9) , J20) . 

Thus, the achievable rates of node s and / are governed by 
this information theoretic generalization of the max-flow-min- 
cut theorem: 

R b s e < mm{R sf ,R M + R c ) 
R c + Rj < Rfo 

Appendix B 

Matching and nonbipartite MWM algorithm 

Consider an undirected graph Q = (V,£) where V denotes 
the set of vertices and £ denotes the set of edges. A matching 
M. is a subset of £ such that e x n e 2 = for e\, e 2 E A4 if 
ei ^ e 2 (23). 

Let x e denote whether an edge e <E E will be selected in 
the matching, i.e., x e can be or 1. Let r e represent the edge 
weights. The maximum weighted matching in a non-bipartite 
graph takes the following form J23"| : 

max r e x e (16a) 

e£E 

S.t. J2 Xe - 1VV G V ( 16b ) 

ee5(u) 

X] < L— J V oddsets U C V (16c) 

eeE(U) 

x e e 6" (16d) 

In ( |16b[ ), <5(t>) denotes the edges connected with node v. In 
( |16c| i, E(U) represents the edges contained in the set U. 
Equation shows that x e takes Boolean values. However, 
Edmonds [7] showed that the we can replace x e € B n by 
x e € 7?.™ and still obtain integral optimal solutions. Thus 
the combinatorial optimization can be converted to a linear 
program. 

Appendix C 
Distributed Local Greedy MWM 

• Each node i knows its adjacent link weights. Node i picks 
the "candidate" node j, based on the heaviest link weight 
and sends an "add" request. 

• Wait for the response from node j. 

• If node i receives an "add" request from node j, i and 
j pick each other as the cooperative pair, i sends "drop" 
request to its other neighbouring nodes. 



• If node i receives a "drop" request from node j, node i 
removes the link from its adjacent edge set. Node 
i goes to the state of step 1. 
The distributed local greedy MWM provides at least 50% 
performance of centralized optimal matching. Distributed local 
greedy MWM requires 0(N 2 ) amount of message passing. 
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